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ABSTRACT 


This thesas analyzes the dynamic stability of submersible vehicles in moons 
in sty degrees of frecdom., A continuation algorithm is ised in onder wo obtain the 
steady state Sélutions in terms of dive plane angle, rudder angie, and longitudinal 
separibon of centers of goivitv/buoyancy. The equabons af mioton are then 
lincanved inthe vicinity at the above stated norunal paint, The eigenvalues of the 
lineanzed system incicare the degree af pabiliry of the worinal motion, Tie 
resulis demonstrate the <uibilicing or destabilizing effects of general three 
dimensional motions a5 opposed to the tailitional tse of straigiit line level flight 
pails. Recommendations for robirst désipn of control Jaws for commanded paths 
in combined honzontal/vertical planes are provided 
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z. TRTRODOCTION 


A. GENERAL 

Az the yarmitiiicy of onderwater unmanned sunmersiples are 
becoming weagined and realized, the question of dynamic 
atéAbpid\ey mist be addressed and justified, Salvage, survey, 
and Tescia OpaTrations 48 well aa detense rajated sperations of 
the future pay become intreasingly dependent on the abi lbicy of 
reliable «ariranne? Undérwirer submersible; co periorm 
prescribed migsions, Theorgeieally, am unmanned submersible 
Vehicle would gor be ae physically restricoed in its motion as 
4 manned yohicle. Secause of 124 dinhabitasts and cargo. a 
monned wubeersible 22 limited moarly Fad arraight line melon. 
This thesis seeno 6o @xplora the dynamic stabilicy mE 
aubmersible vehicles in mecione tn gin degree: of freedom. 

The tight mansovers chat can be performed by an unmanned 
Submersible Vehicle make analysis of oyhamic asabidity very 
aLfibeule. The nonlingaricy of the dynamic propercies iz yery 
comploested and stnecimes very hard fo visualize, The most 
effective analysig sf Ehe Gynamice ond stability of a 
waberbarne yebicle comes from dbserving 4nd/de tiodeliong the 
hyrdrodynamie fortes acting upon i2 Enrough fhe use of che 
Eranplational @nd rotational tquatiens 6£ motion: in six 


degrees of Ereedom. Thes® equacions are linearized an The 


WEeinity of a nominal flight path. The eigenvalue: of che 
dinparized s@t of @guations @re@ on indicator of thé telatiy¥e 
dégeee of dynamic stability for che vehicle's moraion. 

hiven @ mathematical mole: with an established sec of 
physical characteristics, analysis of the Gynamic stabyisty 
can be stedied varying the dive plane angie, Ehe rudder angle, 
dnd che Location Of the center of gravity. Ansliysie of tha 
Linearized aquationg of Sotiog and datermination of che 
eigenvalues by individual talculation would peave tedious. 
Thoretoré,. computer algorithms developed by jJaydal (fet.i| and 
Aydin, \Bef.2| are aged in. this thesis toe accomplieh cheese 
tasks. An explanation of epocoupled stability is provided in 
Chagter IL. Formalation of the linearized equations of merken 
in 3ix Gogredé of fEresdom is presented in Chapeer fft- 
ECnapters II and III sre based and partially paraphrased from 
Papoulias [Ra£.7)]. Tha oxperimental results in Chapter [¥ are 
based on the understanding of Chapters If and IEE. Chapter V 
SUMCizes the results and provides recommendations for turuce 


submersible modelling research. 


B. PARAMETERS 


Severe! Voiriablies were used in the algorithras associated 
with this research. “The shorthand nereation ifs Fram Smith. 


crane and Summey ([Ret.4) 2m is 4% follows: 


x1 Ya @ Distances along the fodty fixed axes. 


ue 7 Translatoonal velocity eaepenenta of made! 
relative ta Fluid along body axes, 


bi ad «& Rotational velocity eamponencs of mode! along 
body acer. 

ei ta ae RyYdrocynam.c force Sempsnents Asong Sody axes 

EK, My ff WaAreky panic mosant compote along Besdy 
oxen, 

yw, 4, @ Yaw, plbch, and roll angles) posicive valor 


Sollowiog the right hand cule 
Ker Fs: Ze hentar of gravity coordinates along beady axes. 


Hees Fes 2D Center of buayancy coordinates along body 
axes | 


Deer pe (De Momenis of inertia About body ase. 


The standard body.faxed, right-hand arthogonal 414 system = 
employed throughout «ll daca simulaciens. All date rug 
pimulatione deseribed herein were inictialilzegd with 4a constant 
propeller speed. of 300 opm, Whach translates té 4 valine of B.0 
Ee/mec for "ui. The data c= analyted for rudder angles 
becween 0 and 20 degreer and stern diving plane angles betwee 
-27 and 20 degrees. Separate data rune wete conducted within 
the above stared parameters for similaced model configurations 
With che center of gravity shifted longitudinally forward and 


aft beereet -1.5 and 1.5 percent of the model Length, 


Lat 


A. HANEVDVERING IN THE HORIZONTAL PLANE 
i. asic Maneuvering 

It 23 S@ktrémei/ impereant fo note thar Ehe Ability tea 
cOntTel @ wemsel unvelyas two Major characteristics ERG are 
Pontrary fo one another. Thame two characreriscics are the 
Sbihity of conreeia Bised course aeabiliey and che ability of 
the vessel fo turn. Obtously 4 Wadeel chat ge exteremncly 
conse stable will not turn #asily, Newton's equations 5f 
motion, used with ceatercence frames both Fixed with syspece ba 
the Barth (inertial) and fixed felacive to the weasel wal 
tllugtrate the dynamics of ship motion. ‘This marian will be 
ftiret described in the horizontal plane and later deseribed in 
the vertical plane. Ae « boty moves within the inercial 


Frac, its Gotions can be described as tol Lowe: 


mi, = Xs (1) 
oP, e F ia] 
t,¥=n (3) 


thtcugh the ese ot the wariables intreduced in Chapter I. The 


SEERESL pe & pPitbains £ce motion tn che itpertial reference 


frame 2 fhe smeoel'’s celativg Pram of reference, Lhune 


CM4T Lions Serome: 


i= ck, coay + ¥ oinw , fa) 
F= ¥, cose - 4, cing . 15) 
x, * UconW - rainy , (ey 
v, = using » v cose , (7) 
and 
eer. (8) 


Therelore, equdtjons jo} and (7) Gam be antegrated and 
DUBAELEWUEOd Inbo equaAglons (1) through (3). Thay ie tier con 
He BUbSCituned inte weecione (44 4nd (2) walle phe following 


Fesulte: 


mi-mer= 2, surge oquatiog, ig} 
ave murs }, seay aguation, tio} 
i, Pre. vate ogakion. (Lz) 


Pest 


These equations 4re trid only 2f the cef@rénce plane center 
ind the Weasela center of Gravity soancide, Cf thé two ore 
indeed located at 4ifferent locatzens, the Pollowing Squations 


Fesatit: 


mM ~ mr - ox rt = X | (12) 
oO = mur + ax = ¥ , da} 
It > mx,(h* url =H, (14) 


poovided that ys = 0) Am impectant potine co mote ia that 
centrifugal forces denote apped4s in equations (4) Ehbengh (3). 
These forces ara not evident until the relative reference 
Blan i Weed ih Squatione (9) through (12), -Cantrituga: 
ferces dre not experienced of the tnettial reference fame 
Sinc® if if sketionscy bit are indeed evident when the 
relative teference plane of the vessel 8 studied, The 
hydrodynamic forces £ and ¥ and che moment W are composed cE 


Eour different types cE forces: 


L, F forogs; Fluid foreés Accing om Ehé uli by che 
subcounding Eluid (water). 
2. A ferees; Contrel surface fotcés due to ridders, bow 


planes, diye planes, thruster, ete. 


i £ forces; Envitonmental forces due ©§6 wind, current 


ey = 
and waves, 


T forces; Propulsion fotees due ta propoallwcs, 


in 
- 


EnPeSters ete. 


ie £6 She CoRseruction characturiatice of ite siterborne 


A= AX, * He tpt T, (15) 
VO Re + tet Fp tLe} 
N©N) + Me * Ne. (17) 


Sinte the hydpedynamic forges Kw YY. amd No dénend on veose! 
Motion 4nd must therefore bé umeticoo of vessel velocity and 


acceleration through ihe saber we can write: 


Xp bal #plu. ¥; OU, iF. £) 4 {igs} 
Yn 7 Flaw, eo. FF i 119) 


Ne * Mp buy We, t 2) (20) 


The incterrelaticoghip cE equacicns {181 chrough |2%) 1s a 
Seplickted ond, Taual maineivering &budiet simplify hrs 
tnterrebationship, Tha Sbgect if fo study Snip resporise about 
a fominal equilibrium point jaqsignated by the subseripe Li _ 
Through che ose cE the Taylor faeries expansion around the 
Riemahal point, and keeping Oniy Che firse order térm, the ¥ 


relationship becomes: 


ay ay, 
Vp = Pelt My Ge Met hy) + fam By) a= eiv= Ml 


ayy (a2) 
“iF 


+ ee (f- FI sS 


where all of the partial derivatives ore avaluated Ar 4 
Sominal set of ‘conditions. x, and , can be Fimiisariy 
minipuilated 4nd avaluaced., 

The sominal eat of canditaons is the set of variables 


associated with straight line moolon: 


i il ee ee? ie (22) 


For tha ake of inderstanding, the vesdel t= considered te be 
Part iatarbaard aymmeiric. Therefore OF / dy «© OFA) «© 0 since 


4 change ino the fosward vedtocity os acceleration will produze 
no transverad Forde in Vessels Chal are salad si¥reecric in x= 
plane If goditden, 22 the vessel le indeed ip @quilisrium, 
and in Ar&edaghe ithe metic. cthiftte can be of throes. 
Therefore the firat CSrm in. @quaticn l2ii is 6i40 sara, One 
these additaondél Considerations are, made, equation IZ) 


PeWuces bos 


In J2k® Saghiad, thé surge forcé aid tie yaw moment can be 


SIMplLitad: 


ax ax, ay, ay, 
= Fie ak - = F + F co - = 
Rr? 3a." rth Oe” GEE a aF¢? 184) 
ay, an ay day 
Ny, * = ¥ - ag * ad + art : (25h 


Where the cross coupling decivacives 4Y,/d°. AY, /08: Any fd, 
afd Ano Wouelly heave amall nonzero Valie® betguge mage 
VESSELS are ot symmetrical about The yt prent @vee we thet 
plane is at the midiengto oF che veesel (‘bow and stern 


CORtOUrS are Usually Grastircally ditferesci. Hiwevyer the 


Croas coupling derivatives &ifdr, 4%,/45,. d¢,/dv, and 44,/aF 


are =erc because of the ayreecry about the «2 plase and 


‘eonbtson (2a). EodAtion '24] sow reduces bo: 
ax, ay, 
R= —fi + Slips ; i2Gh 
eu u tas By 
Using stepdard pobaticn., fide = Vi. GHear =< BH. ..- SEE: 


Thago are known as the hydrodynamic derivatives, “The physical 
meaning oF thidd mathemacical dertvatives ia very imporbant bo 
the deyelopment of ndérstanding the confection Serween 
maneuvering and hydrodynamic forces, For axvample, ¥, ie che 
force in the sway tirection due cto a unit change ir che yaw 
angular Wolocity. 

Using €has conventlonal nobarion, and subeclturing into 
G@quations (13), (13), and (14). the Jinger equations: of pester 
in the horizontal plane in the absence of environmental 
disticbances (idd@aél handling conditzens! and with the control 


aurtaces fet af zerd, become; 


(a= Fd *elu- oo, (27) 


(m= Fj)o- {(¥,-mxlt= Fee (Ye amr (28) 


(J, - 82 = (X, - mele ae CM, - oir. (25) 


Wire we Gan see that Witnin Lineairity the Surge BqU4ATion 
decouples From away find \aW, 

Hote that all of the cerns of equanionm (28) and 129) soe 
ing lude the effece caused by holding che vessel's rudder at 
zero in orde; for Ass tht previously; itimpostd. simplifying 
conditions fo hold true, = che Pomsibrlicy of now-contrals 
timed mocion were fo be considered, equations (28) and |29) 
Aust include terms on Cheie TOspeccive Tight sides expressing 
the fortes and mement= caused by control 4urtaces déflect con 
(i,@,. Eddder defleceson) 49 4 Function of time. Assuming thar 
the control surface sear generally used on the horizopcal 
plane io indoed « tutor and that the sadder force and memenc 
30 the vessel 4ré@ functions of the ocodder angie & eniv, ane 


noc A, tha fFallowing §s Erua: 


(0) = 4, (=o) + al =O (30) 
ay, . 
Hp(B) ~ ¥(8 Oh = 28 = HS, (31) 
aN, 


WiylB) = NiO = 0) = SS + m8 (32) 


Pogibiye tidder deflection carresposd= cto 4 Fart OPsrm when 
agsociaced with «a @Gerh rudder YY, and We are ‘the cudcer 
hydrodynamic derivatives, Tf che rudder forces ond pomeocs 
426 trow 2fcliided if the dinwdaricec eaquactons for away anc yaw 


jgquations (28) and (29), respéctivelyri, the equations of 


BEtian become: 


im-¥Jt- (y,-m,fe yee iy, - atjre ye | ON) 


Uf = Nt = (My) — mor) = oe = UN, — meine + yd 0) 


3. Stability of Straight Line Motion 

The concept of an abject remaining of & g3@#f path is 
related £o the concept of 6 \ae8el's course grabilicy or the 
snabidircy oF its Straight line morion. There are varios 
Kynes of motion stability associated with marine vehicles and 
Ehey dre separated by ehe aterifurces of their ifitiadl 
@plilibraem state that are retained in che final pach of cheir 
center of pravibys These types of motion atability are 
Gepicted in Figuré i. For eich type of stability depicted, 
the Yeesel is initially travelling af constant speed alang «4 
Straight path. Case T, straight [ing a= dynamie stabi licy, 


retains che attribute of straight tine metipo but Abt the 


—_—_=_— —_— — 


ne 


MTEMOED Pate 


So TAL fare 


He 
> ier oF SoTuAnance 


Lfuire 


Various Sipds of Medien stability 


[ASE 
Ta Sed FS kL 


Sah aim 
eesti IF sein Ty 
Lo ars 


Sake ae 
ecm osa). tke 
ape AAG 


Saas on 
Se TIOvAL TTAMILIT® 


athripute of Direction From the iniaiiai sate of @quilibr: 
4tter release from a disturbance. In Osse If. directional 
Stebriity, che fentl path afrer teleaue from os disturbance 
Tetalne nat only the sezaignt itne attribute OF Eho initsal 
path, | but alee its dimeetion, Cage Iza is similar €a Smee IL 
except that the vessel tops not ogectilos@ otter fhe 
distucbance, buck Passes amoorhiy to the same path AS in Case 
tI, Finally, in Case IV, positional motzon gtability, the 
vessel returns to the original path ant retains direcrcion as 
Wall 66 POFLtian, 

The typee af atabilicy are clasaified in an ascending 
agiger, A Gitectionslly stable vessel mist also podads 
Straight line atability. One that 2r@ positionally stable mist 
posses both directional and straight ling stability, Straighr 
iine stabilicy or instability is indicaced by che solution to 
& second ofdér differential equation, dirécraonal stabiliry on 
inarabiliry ig indicated by che solution vo 6 third order 
differential equation, and poritional stability or instability 
is indicated by che aclutiean to a fourth order differential 
Equation. 

Each type of controle fixed acabiliry hag am agaociatead 
numerical index whieh by itt gign detecsings che felacive 
stabiliry ot dnostability oft che ship: The magnitude of cha 
index determines che degree of stability of apetabiliry, im 
order co illustrate the determination of thig index, and to be 


able to evaliate the degree of dynamic stabiidicy, a rebar bo 


i+ 


ithe cdiftierepetal #aqudtiene of mition f28| and fo) 48 
nétdeaary., Thig bomegeneoos get of aifferintial #quarione 
Will Mave @ solucion of the form: 
vich = ya * ea (35) 
rte) = ce © ret (36) 


SHEESH CHE SOAS EARES Vy, Mya Ty, Fer are determined from the 


boitaial TIndi fc ioensg . OLEfareanciacion ang Eubsequent 
WuUnTEIEURIen BE equations (94) amd (a7) inte equationg (26 


HAG (29) yislos the following characteristic squarion; 


Avie Bp = =o , {37} 


wih 


A* 12, - Mp) lie= 4 - (¥, ~ meg) Oy - mel; (38) 


Be =U, =.) ¥,.- lm - ¥) Ue - met 


= (¥, - mt (Ny - meh = UF, = mee (35) 
c= (a; - mxglly, — (Y, - aulM, , (40) 

Eth. Tages 
ayo, > 31-2 2 (Bi 01S) can) 


Th Gn ciurse stability is fo be diaplayed, then co. and o& aill 
have fegative cmal parts. Seferring back to wquatiens 176) 
and [327)) Bork viel and rle) will go te tere for tocreasiog 
Enz . Tats hand that the ship will mvantualiy eassiime a4 
Straight line direction at a differant, in general heading, 
This j= straight Line stability as defined in Figure i, 


Symbalicaily, the stability requirement 


Rig,, a) = Da, 


tranaiates te 


a) 


with all ehe coeftiacitente of bhe quadratic equation |MEl are 
posiciye, To evaluate fhasée condiciones, wt Gilet <héek Ea 


signe Of che hydrodynamic dorlyatiyas that fice up a,8 4nd | 


A physica d iptercprecation of che hyvdradynami® vattianlie= as 


force: and moments acting on the hull dow co the varicus hist 
Motions with respect to che water is necessary. 


ALber some andlyeig, the following conclusions are drawn: 


Always negative and larga: 
always negative and large, 
Always negative and large, 
jlways negative and iarge, 
Small and 4! uncertain ssgn. 
emill and of Wntertiin sign, 
emali and WHESrEAiD soon, 
small and Gi uncertain sign, 
stall and of uincertagn aign, 
always pPosLtive and large, 


always positive dnd Large. 


With Ehese conditions described thusly, A and 4 sr therefore 


Q1lwivs BPOSth ive and stabililzy must require: 


ce UN, © mxgliy, - (FY, —aitiW,* Oo . (42) 
This is Em Eeabilisy Criterion, The tore positive C a2, che 
more stable the ship will bereme ond the mere increasingly 
diffieuic the ship will be ta maneuver, The more negative ¢ 
=¥, ih@ wore uUnstabla the shig will Sécome and the cantinudus 
vee 8! rudder will be required to maintain course, &ltheugh 
% 15 aeall and of uncertain sign, for stabtliey it ts desiced 
bo have 3 = 7 ieetitér oF gravity Forward of setdshipasi. 


Stamdard shir design #xperience agrees With this ohservarion, 


H. MANEUVERING IN THE VERTICAL PLANE 

1, QDdinear Equations in the Vertical Flane 

HaneuVers in Che Vertical plane can be deacribe in mich 
the S47 manner 42s they can if the horizontal plane. The 
degical pregresaion in the dererminacion of arabilicy wall be 
Emaced out here for the vertical plane as if was ih the 
prévioua g@ctlon for the horizontal plane. 

a Submerged veggel moving in the vertical plane has 4 
nominal path defined ae the Level pach 4t the commanded dapte. 
Tt has «a pitch angle, 8, whach is defined 4s the angle berween 
the Vehicle's longitudindl axis and the x aki and an angie ot 
aceack, @, which if dwfined ac the angle between the » 4xia 
and the total velocity vector, Figure 2 dllostrates these 


POIngs, Ii Newhon*& law of Botion 9 once Agaan uaed, the 


Le 


i 


Figure 2: Vebiclé Geometry in the Vertical Plane 


SSLLOWLAG SQUAELOLE S22 flue Using Che variables: Persducey is 


Chapter I, 


a i 
mi; = Ze 4 


sh =a , 


(a3) 


Pai) 


1e8) 


Af 20 Tne Mor Sontal plane, hs vertieal trankéermaction of 


SOoRGLAArPEe System Yiehag: 


x, = +kcos + taing 
2, "= ~Xsin® = Zcona . 145} 
7 sho = wsdnh | 


2, = ~usinh + smosi 


whese the equéteans Of meckon ape: 


n-ne = £ , surge equacion , 


m= tug £ , heave equacion , (a7) 


ia" , pitch equation 


Ino 


feg. ia} 


When EYCresped with respecciem the vesgel's selerence trom 


these equeticug becca; 


my my - mk + oy ey. (a3) 


26 


we My - ey - og = (50) 
yg = @ix,( = ug) ‘mx li- wg) =, (51h 
Wwheee (%,,4)) re ethe coobditiates of the ship's center ct 
qraviry with swapece 66 the ship [ised reference frame, and we 
Assume i, = U. 

Equations (49) thraugh (51) demonerrace 4 coupling perweer 
Surge sorton and verbicel planes rotion that wad not preaent io 
the dipcuseron on the horizontal plane, TE the vertical plane 
M@quAE Lone o£ motion are linearized around a Level Flight path 
a: gonstast speed, U, we can see from equations (45) and |51) 
thac th® wh berm mH 4nd mn will teamain afte: 
Linearieation,. This means that heave and pitch are coupled to 
surge #ven far gmail pecicns. This dynamic coupling arisea 
due to the ooocero = berm. Whe a vessel pitches, 3a nonzero 
gurge @cceleration will contribite Eo Ehe tnagrtial momdnt. 
Shyioualy Chere wae no such linearized dynamic coupling in the 
Hormzzantal plane. Altheugh vertbeal plane motion must br 
Studied cogecher with surge, the coopling #iscusaed here is 
cabwtively, gmail due to small-yaluesa of and in Enis case 
ean be neglected, Theretiorg, the uncoupled Linearized 


gquatiane of mation in the vertical plane are: 


aw mii ~ mya = =. (S24 


re 
os 


fa — mkgtit = Ugh sae: (53) 
hich Like in che previgua Séeclon, = ond M cac be expressed in 


terms cf slow amecken dirivatives., Vaitg vertical plane 


adtation they ara: 


Ee Sw + fyr+ Zw * yp + 28. (Sa) 


N= Mowe Mags Mw Map fe, - 2) oinb «4d, (55) 


where Z, stdnda bor the change in the hydredynamie force in Z 
WLEH @ Unit change in w jnotei Z, £8 synonymous to ¥, in the 
ROPezontas Blake! 42nd tho sao for che other variabies. ‘The 
additional term in equation i585) 16 che hydrodynamic rescoring 
moment, @rnd ol) other veriabier are in accordance with the 
ligting given on Chaprer fT. Peoavided the vessel is properly 
balanced (weight «= buoyancy!) Ehite io no auch hydrostatic 
force in She heave omiation (54). Tha & im Chese equations 
Sefera £5 the aterm plane angular deflwetion 4s shown in 
rigure 2. 

Equatiota (52) chirsugh (35) iéad ta the .inesar uncoupled 


equationg of Motion in che yertical piane- 


lin Z,Jw— (mgs ZG" 20 * (f+ mig = eh, | 156 


2a 


(2, * Mgh@~ (mx, + Mw = M+ UM - mtg — (2g > Zyl + 8 
(57) 


Thees Ewo equetians can be ised with G@quation (48) LA otdear ro 
(iredier the response oF Che Veasel an El@® Micinity oF a 
nominal devel Elaight pack. Deviations in oceprh Bain be 


Sampucec Proms 


fF © -Gein® « gcosd , (SB) 


or 


#£e=78 (59) 


‘8 des linedtized form: 
2. Stability of Hotion 
SEGbility GF motion in fhe Vertical plana can Fé 
atudied by using equations (56) and (7) after & is gear te 


2669. Ag in the cage of <he Gorteonta: plane, small marian 


4 


are being Sthdlec Arsund the nenineal Fiight gaith wich Ebe 


eontcols Fixed. Equacsshs (34) spd ($7) beeome: 


fo - 2,)8- (oxo Z.)e- owed Bg. (60) 


(I, - Myber-imx, + Mw Mw tM xg = yee oo, C82 


where 


ty = = (ty ? Z,/W= wits , zy) mag e (ez) 


and = 1 cine: dimensianless units are coreoly uoed,. Th) 
tint 24 Ele hydreskaAbic “Spring constant* dud te Ehe poqzera 
MeEACenEfic height. Tha <orresponding term in the horizontal 
plane be MH, and be eqdal £6 Zero: The negative sign in 
equacion) (63) cakes MH, = 0, which is consistent with 
@steablished definitions for positive forces and moment: in the 
Fight handed orthogonal toordinate syatem, Manipulating 
mqugcions i601, led), eamk (06) snd ecransfarning ta che = 


domain gives the £1 lowing: 


(a= Ze = we (=e, ema = (2° male, 16a) 


|= (N+ mela - Ml weel it) = wa -(N, - oxe - QIE =, 
(oa) 


where he charactoristic #queationg ga: 


Ag'= Agi +f =D=eo , iS) 


whire 


A= la— Kill, - My) - (My * Gayl iz, * mx) (66) 


Be -(m — 2) (i, — mtg) — (ty = mx), — (2, + mbt. = mx) 


“(dy — My, (67) 
ce (Wy - mle, - (2,4 mM, - de - Mm (68) 
B= 2m (69) 


When coMmparéiag equétions (86) through (29) with equations | )H) 
Ehraugh do), @oefficients A thraugh C @re simular itn the 
horagontal ard «wertical plane@a With the @xtra matacentrie 
{Hone Term Appearicog in equations /65) anil /afi: Ib there 
were oo hydrpstatic teres, thts cubic characteristac equation 
f65) would raduce o6 @ quadratic similar to che equation 
derived 35 the tericontad plane, THis means thar there 48 « 


Earo peat in Ehas case. A Piysickl analyeis 25 neeasAry to 


understate the bul!t meanieg of thia mathematical face, In the 


Vertical planes 


{Ag! + ds! + og + OB = 0 (70) 


This means that, Aadliming oortain conditions of stabslity are 
Sacistied, it i# possible for 8lt) +4 ae cs, Therefore, 
Gireccional arabalicy ts posible. Where in Ene herisencal 


zisnc, 


(ast + Bo * Chap = (de? + Be OP = 0 ttl 


This meant that Ene fest Fhat can. Be expected in che 
horizontal plane is Piel = rie) +0 go be, Which moana thee 
w wild be Gnergasing linearly with time. Therefore, 
ditectional stability is mac possible and che best possible 
SILMBE20n. 2% strenght tine stabilicy. 

in the cage of the cubic equation the Rourh-Huricz 
condiczons for srabilicy dictabe that all cosificients A, B, 
C. and 0 be positive and Ehar: 


BC-AD? Oo .. (Ta) 
AfRer close analysis simittar fa Est sanducted th the 
Hearcizopeal plane, “Ut aan be proven that A. 2 and D «are 
fcaicive, Thereforé im order for che Vessel bo stable in Ene 


vertical plane: 


C= (M,- mez, -(Z, + mM, -im= tm) 0. (73) 


Slace jm = 2) = 7 and M& = GO, this comditien ia saeler toa 
aatisfy in fhe vertical Plane chah in the herisental plane. 
This proves feat ¢ positive mecaceqcric height cantributes ta 
aynamic stabiiity as wall os static stabilicy, Assuming that 
elation (7]) is satiafie\ed, if cain be shown that equation (TZ) 
i5 Always Satisfied, This is best shown by using a sechod 


derailed ty [Ref.3). Tf the coefticianes are written as: 


where the subscript U retersa to the bydrodynamic coctlicients 


without the Hydrostatic term, and subscript 1 sefercs to the 


MECACRREPLC twee. The pirameter & ce the metacentric Aenzghe 


te Zl: 0 Therefore: 


BC - AD = BC, + e(8,0, - A.D). (74) 


Based on the given conditions, &, ,. A. and DB, are Ali 
pasitiva. Therefiote tha difference |Bic, Agly} is fH amall 
number Of uncertain Sign. Simm it be multiplied by c, 4 
amaui positave mumiber, the cesalt af ao small «umber ar 
uncéernain sign. This small numper, regardless of sig wilh not 
change che figd if @qlsation (74d), assuming that is 
positive. Tiis means char © ©» GO will be «a sufficient 
condition te sabisiy vertical plane stability. The vertical 


plane 4cabilipy peducea ta: 


(mM, ~ mx) Z, lz, * mM, 0 (75) 


Whore Ehis i¢ 6 SUEFiCient condition for vertical atability, 
4@ if oc 4 feressary condition. Tt is posaible chac the 
dradynamie becom will stabilities the veooel and equation (72) 
will be satitsiued fand all coefficienr= will be positive, 
even though equacior (75) sight nor hold true. This. condition 


La dependent on the magnitude of che metacentric height and 


che apeed of the vesse: ha & Veseel thodelencos an incredae 
in forward apudd, 1t Wilk behave ipsa atdoly in the Wertieal 
pLane, 

Aessiieng the frabilivy Criterion tos the berber 
pidnn Etc @guAtion (i42), equation (|75) san sé2ve 4g the 


Stabiltey Crike 


Py 
E 
a 
3 


for che vertical plane. These Ne 


AXPTRESLORS are Commonly nhormaliced ina Stabilrey igi cours 


gpk, = lg (76) 


EL(M, -m,) * 


NotY, = mh 


2 ———— iv} 
¥,UM, = mors} 


These Ewa Lndiees 4re thé weriical and horigencal plane 
BEEDLLICY LAGLeAn ane They off Used fo asses feneliverabiiicy 
and selative srability. “They dfe Slways less than ona, and 
Boost ate ‘alums indicste stability and segarive values 


indicate pistability, 


2G 


Zit. MWANEDVERTHG It SIX DEGREES OF FREEDOM 


RIGID BODY EQUATIONS OF HOTION aHD KINEMATICS 


Pia 
BeveLoomence of Ene general min@uvering eguationg 25 Six 


Gagrees aoc freeocem requires Che employment of two Geeraisaéte 


systems depictet tn Figure 3: 


Sik 4.81 | twEerente =herkial Stame, 
a ae ee moving, feare fitted sn cha vessel 4 


SOonel ssf =Shter , 


a Ty 


nd gody 


=<" 


Figure 2: Fismd anh Moving Csotdinace diest dee 


an 


As discussed ip the previous chapter, Hewtoo's law requires an 
ipertial frome, whele the m@motdon of the Weasel ya best 
destcrpoed in its own frame, & dendtes the center oF gravity 
ef the moring tinny, whith £8 ipbated af (&,.4,,2,) with reaspece 
fo Che Alfk,y, =!) Eraew®, Unat vecborp in. this moving frame are 
denoted by f.9,38, an the x.y. Glrectians,. respectively. 

Tha translational equations of motion Can te faund chrsugh 
the ute of the principle of Jinear momencum im che inerctal 
trame whe pate of change of the Jinwar maqentum equals the 
Sim GF the external forces. Vector dnalyele le than ised ta 
expreéy these guantities tn the Av eystem, Chapter 4,1 of 
Haference ] explaing Enis conversion §n decabl_ ‘ace Epa 
ronviréoon bo ERe@ A siebem i= Accomplished, Mewton!'s law car 


be Broken ints Componanrcl Gat 


m(ad + gw - re - asia? 4 3") + yi ipg- ft) <a ipr=@)] =x , 


mle + cu —pw-yalr* =p") + aiar- 2) = alaps fH) =F . 


af = pv - gu-a(p4« gi) -«aicp- a) eyolrase Bi =z , 


wiere Che Soka: exciting force .F, ts denored yi 


Pext- f= 2k. (aL) 


THereiore fhe SXternAd moment. can be wriccen as: 


Mext-~ePiwk , rezi 


“He ALOkee LAS SoeponGnls Ags 


ig6° tf,—1ler + Ler g - tte - 24 - tp = #} 
-m(¥glie~ug~ vp) - ele ur=wp)l «=e. (83) 

Rags ft, ~ ior - ft tge = pl + tofag- 2) efile — 2 
“ala lw - ug vol --agld - vroe wl] = A. (oa) 


S,f = (2) - RIB igi = 2) - 7, ipo fd) + Tg = 


tn(xs( + or = ep) = Pgll = ve = ey] ev, (85) 


Equation«® (74) through (40) are the translational equations of 
motion, @gnd segethar with the cocational equacions (43) 
through 85), they make Up the 2ix degreéa OF Ereedom of 


motion of 4 rigid Bady expressed in a coordinate eyatem moving 


42 


Slong with the body. oOntt 4Sg8in Ehe Wariabled involved are 
the sage Veiriables intredudmd in Chapber i. 

Although the equaticgas of moitias are weitten in « Erame 
Associated with che vwerehicle. the veansel’s ortentacbon and 
poeatdon bn the inertial frame ia atlll avery important. The 
AnQULAT (OrLentatuagqal Telationship Setweoq *hoa two Eramee ia 


determined by the three Euler angles 


yaw of heading. w 4 Fotation Sbout x exis 
piteh of ealevaticen,& | THOtatgon about y a2#15 
roll of bank, @ 7 TotAtion AbOUE = axis 


Amy vector dogceelbed in bhe feo sceferance frames ora related 
by Some Trans ECTEACION matrix, T which ta a funerton of Eheae 


throc Euler angles. 


fnaceiare, 


re T(¢.6.¢) . (86) 


Tha geanoreal eranstormacion marrisx can bw oroken up inte Eres 
individual robabioog. The order of application is very 


Eribtical in offer §5 maintain the jmbegrity of che 


Lranslortmation. Sny@ically, <sniinite cotatiens -de cof. 
Conte. The prescribad manner i# to first aw. then pioch, 
and Finally t& sell. Sbandacd optriv poration therafere 


ya Gi Tia) 


THC,O,g) = Tid TIBI Tig? . (874 


The Gull fransferssriow Mibsik i=: 


cosfconw cosfsiny -sinb 
Tid, 8,4) =|-sinbeoabesingsinfeosy cosécosbrsindainbsin singcosh | (68) 
singsing+cos$cosysind -sinécosp+cospsinOsiny cosdcost 


— a 


Sal. Perse ions do commute, and tor mall angles, 


ca | 
Tid,8,¢) =) 2 @) - (63) 
6 -# 1 


ond in this case Ti@.8.@) = Tih, 8.@). Hquationa (47) and (46) 


Ch0 few be used Co orang fore Amy Vector if che them ba tk 


44. 


aqguiyvAlene Ly, Che Gener fatem Chirough Eh lee <r makrex 
Bultrplicars sn. Here 4fre Some ieportant feampies AF fhe 


relactonshap between Eno Ewo teference Esamep: 


#, ® uy cosfeos + ¥i-sindensd © sindairiroa¢) 
+w{ainfeiod + cosbeoswetna) 30) 


Y= UcosGeing + vicosgcos¢ + sindsinbsiny!) 


‘Hi =-Bingcong © cospeinining) , 191) 
5, = -u dinl + ¥ aingcos8 = » coséconl (92) 
pr -psind «6 , (93) 
g= 4 sindeosh +8 cond , (94) 
re cosfcosp - 8 sind , (95) 
$= p+ geinbtan8 + rceos$rand , (96) 
Begeosb-feind , (a7) 
; sing . ,co 
VY Doe eae aoe 
Mac ce Sart oo FePPerenkS fhe .netrass. =oPf erence 


frame. Por amit angles, cquantons (¢4| ehooudn (98) become: 


ep. beg, beer, 


which meas thac the ratea of qhange of the Euler angles are 
che same af the 4nguias welecities only for infinitesimal 
rotations. The six Kinemacsc equations |90) threugh )92) ana 
196) Efreugh |98) aleng with equations |/7H!| through (80) ang 
(44) through (#5) provide a coopletce descripeccon of che Sod 's 


HHCLSo iff six tegrees of Iryadam. 


BH. EXTERHAL FORCES AND HOMENTS 

Equations (78| threugh (80) 4nd (82) through (82) are the 
Six gégréées of Ereedom equations of motion for a ragad Body 
wherow the left hand terms: represent inertial forces and 
mamnencte while che right hand terms represenc che eaxternal 
forces And moments: “head external forces and moments ifclude 
hydrodynamic Scontrsbutiaons ge well 65 weight snd Suoyancy 
effects and forces from contral surface detleceions and 
Propeli@r Ehrust. tydrodynimic festering forces and mamenes 
Sree due Fo FR@ iehitle weighE and fudyancy. Forces and 
TOES Gu To contri Surface deflections are experienced as 
added drag cn surge, while ln sway, Heave, pitch, and yaw chey 
Are directly proportzonal to the deflectien: dust 4s he 


position and valocity yectors in the previous section oan be 
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Gipresaaed in terms of wither the inertiai reference frame or 
tha vesael!a relative reference frame, a5 can cheese Forces and 
namence, She cre betwen Efe two Prams of roEwrence in, once 
again Fhe toandformaclon sacrix, Particularly in the case of 
fhe hydrodynamic Eoraf and jonent <onteinutions, the 
derivation of thease binding expressions can be very 
complicated ard lengehy. :£ a morW in depth discussion of che 
subject if deeired, Reference ) prevides an excellent 


overview, 


©. LIBEANTZATION 

LEfGAFLZation, @& EaUGHE Oe & mathemitical pad] 15 Useful 
for untierstanding the Behavior of nonlingar aystens. The 
principle ia used here to analyze particular vehbele nonlinear 
mation 4bout & nominal Plight path. Generally, Linoariastion 
16 performad in Che vieinity of @ somlneal point. The poine 
is that point where the system is expected to spend mest of 
ita life. Te Ootinal point if) defined by che =yaren 
condition that ail time derivatives equal sere. Any existing 
syetem contrals can be faxed Gr oxpresecii ag 4 EunceLaon ef the 
fominal point: Gsce a syatem aolut ion has bee computed using 
the nominal point and the necessary qualiftiog condition, che 
System can be jinearazed. Thie in done by determining the 
Jacobian mathices GE the system. This sew linearized system 
can be used to anAbyé syater apes loop stability thropgh the 


axamingtian o£ tha sigenyalugs of the dacobian, 


IV. DISCUSSION AND RESULTS 


A. GENERAL 

The major motivation GE this chesis ig ce firther explore 
the bine Gi raasoning that underwater yehicleas are ao Lenger 
bound Sy the atracght line mobien chat is caused by personne! 
or cargo staged onboard phe weasel, This thesia sets tO 
pravefverify that 4 particular underwater vehicle tis 
stable/unstakle when enough physical lofermation 18 Snow. 
Wizie gathéring infermation co form the stabLiLtty aAnalyo:s. 
Yaluable informacion can be obhainesd concerning EAe morien of 
Ehe vohicle while Gt performs three Gimenticnal maneuvers: 
The objective of chis study 18 Eo provide 4 paramecric 
ghabilicy analysis of the Swimmer Delivery Vehiole in six 
degres SE Eréetom motion. ALL phiyisleal prépertiee of tha 
Wohiela 628 ‘now «¢ well oe che [inssar, Hanlinesr and 
bbupling terms. The main Sbhject i9 co aases ooncbol design 
while Keeping the controls. Sixed and varying the distance 
between che longitudinal center «ft gravity (LoG! and tha 


Langitudinal 6£& busyaney (LEB). 


BE. PROCEDURE 
AFlir & Drerequasite understanding of uncoupled motion ith 
che Herizesniteal ond vertical plans, the general procedure 


followed in chis paramecric study continued with Ehe full six 


aa 


degrea of Ereedom equaticns st mobic, The neteadacy 
variablec in order to andlyse Efe system oynamice ace the 
‘rindglationalL velecitiag surge, =way, afd heave |u, WW, wi] 
Tetetional velocities roll, pitch, and yaw if. G) lp the 
Euler angles OF roll and pitch (6, Gls and the cate af change 
f heading - the first derivative of the yaw angle. @, in 
@rder ta provide a comprehentive analysis of vehicle sehaviar 
and capahilicy, a wide vatgety of realistic control parameters 
and loading Gonditiong were ¢hoséen and evaluated for thaiz 
GOneEribytion fo che steady state morions of the above atared 
Wairtablet. The final scep in che proegdure ie to Iinearice in 
the viciniby of the wteady arate and to compute the system 
eigenvaitien. Thtse €19@nVdlun— Shryo ds 2andicators of the 


rejacive stability of the yahicla under che given copdiftiona, 


C. STEADY STATE ANALYSIS 

The jsateady stare analysis Weerted with che Fuel, 
Honlinear, coupled six degras of Fretdom equations of morion. 
All cise gerivatives were get co 2ere since gsepady state i= 
detined ag thet point whare Chere te no chaonge in Che orate af 
a waT?éble With réapece fo Eis@. The resulting 4lgebragec 
ayebem GF ¢oaupled, ponlipgar myuation# was solved for the 
conbroal paréuatera and loading conditions tsed. In cider ca 
pectorm am adequate analels,, 2 Cantinudtion slgorishm 


developed by Seylel (Refi!) , with gubroutanes bop Aydin (feF.c| 


werd waed; ‘The Tellowing sange of contcol paransters were 
asAnineds: 


- Greco dies plage angle varied herween =) 
degrked 

= fuader angle varied between 7 and 28 dagress 

- LOG/LCR separation berween 21.5 percenc of vehicle 
length 


The graphical resulta’ Aré opganized in tets accétding tc 
LEG/LOB separation. Thé LOG/LOE separation Was varied between 
41-5 -partent 6f wehicle langthein increments of 0.35 percent - 
Each fer CE the toliowing centinuatian results depioct® the 
Sntire gtere diving Glare angle range as the = axis. Each s3e@F 
has 46. plot where Ehe y Skis is oné of Ehe Afcrementioned 
vVarigbles go, ¥, my By gy. =. 8 8, and F, The translational 
velocities are plotted an fi/ses, and the angular velocreles 
in gad/aec. ALL similations were run at 4 constant propeller 
spand of S00 tpm, which translates to a value of 6,0 €t/sec 
for the forward apeed, U- The final plot in sach set is the 
telative degree oF srabilicy whieh wilt be Hiseuseed in che 
HEKE S@ctien, Theae plots mist. be Analysed as a Set in order 
to provide @ compreliensive lnderscanding of the vehicle six 
degree ct freedam m@tien. Throwgsh the somparisan of these 
plots, thé metion Associated with asy maneuver within the 


range of parameters can be predicced. 


Ah 


D. STABILITY ANALYSIS 

Steady state analysis GE the manguver alone if ooOcC 
suttgcpeant, Lt mist be determined :f thea yihicle La to rien 
stable (controllable! whiae cha maineuwyir |5 fegoriacnd. 
Therefore the final ploc of che sec depiers the degree of 
Shabilipy plotted aghinse the range of the stem drying pias 
angle. The units oE che degree of stability ar@ un terms of 
invyerso t2ne, gec'. This index be obtained by evaluating the 
4 4 9 Jacobian matrix opeated by che determination Gt the 
cheady atetes diacdased in the previews S6ctian, The 2argest 
rasl pact of all resulting eigenyslues an &@ cCentinuation ia 
che degres of stability, This, in physica! terms, |= the most 
dominance time gonstant tor the system. The indies measures the 
GigeesE Exponential canyvergence fo the steady staé&te sation 
when negative of the Eastest axponentiad divergence fram che 
Steady 4atare motion when poeitive. The S58 wegative the 
index, the more ateablée the Vehicle mantuve:., Conversely, the 


fore® POSLELve Ehe Lidex, che more unerable the padeuyer. 


=, GRAPHICAL RESULTS 

The following figirés 4long with an understanding of cine 
previoclia sections of thia chapter tilustrate to the roader the 
aynamic# 6: yehicle motion and the gfabiliey of that motion 
within Eh@ paraowters inputted into che fontinddarion 


iigerithm, 
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Figure 7: po Va Stern Plane Angie: x= 08 
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Figure il: O vs Stern Slane Angle: «y= 0% 
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Figure 43; Stabilicy Index va Stern Plane Angle: == i 
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Figure 21: 8 va Stach Flane angle; wy«.34 
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Figure 23; Sta0iL2ey¥ index va Stern Flane Angle: «8.5% 
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Figure 24: u VS Stecn Plane angle; xa=1.0W 
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Figure 25: v ve Stern Blane angle: x yzsl.d% 
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Figure 27; 6 va Stecn Plane Angle; x. 21.08 
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Figure 28: g vs Stecn Flané Angie: w,=1 0% 
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Figuré 31: @ vs Stern Plane angle: x,=1.0% 
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Figure 32: Gf vs Stern Fland angle: =.=) .0% 
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Figura 33; Scabilicy Index ve Stern Plane Angle: xy=l.08 
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Pigure 35: v ¥5 Stern Flané Angla; %y=1.54% 
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Figure 37: p va 2temn Flan# Angle: «41.54% 
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Figure 323: 9 vo Stach flane Angle; x=1.54 
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Figure d1; 9 we Stern Plane Angle: 4.41.54 
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Figure 45: vovs Stern Plane Angle; x4=-.3% 
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Figura €6: wows Stern Plane angle; x{=-.5% 
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Figure 47: 6 ¥H Zeer Flane Angle: x.2-.28 
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Figure 48: gq va Stern Plane Angle; #,=-.54 
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Figure 49: = vs Stern Plane Angle: *s-.34 
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Figura 50; 0 V3 Stern Plane angle; «.=+.5% 
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Figure 53: Stabilicy Index vs Stern Plane Angle: i,.3-.54 
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Figure 54: U'vs Stern Plane Angle: s,s-l.U8 


Figure 55: v ¥a Stern Plane Angig; Sye-1.0% 
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Figure 56: wows Stern Flane Angle: wy=-1,0% 
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Figure §7: povs Stern Plane Angle; ,=s-1.0t 
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Figura 59: ££ ve Stern Plane angle: +ya-1 7% 
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Figure 60; @ va Stern Flane Angle: «,@-1.048 


Figure 61; 8 vs Stern Plane angle; #j,)=-1.94 
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Figure 63; 5Stap.lits 
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Figure 64: u va Stern Plane Angle; x,="i. 34 


Figure 65: v vs Stern Flane Angle: “,_e-1.54 
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Figure 66: wove Star Flane Angle; HX g<-1,34 
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Figure 67: 7 vs Stecn Plane Anglia; i 
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Figure 68: q va Stern Flane Angie: «y=-1 55 
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Figure 70: 6 va Stern Plane Angle: «,=-1.54 
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Figure Tl: @ va Stern Plane Angle; =,=-1.5% 
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Figuss 73; atabilicy Index ve Stern Plane Angla; x,,=-1.5% 
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¥. CORCLOUSIONS AND BRECOMMEDATIONS 


The results GE the previous chapter show that the 
COnTLBUetion Algerathm 2a efficent in eracing Che ofeplical en 
SLeady SCALe solUtLonk dag0e1ated Sith mations 1f sik degrees 
ef Erecdon. The iiranslacional and angular vedlocizy 
eanponents; the Euler angle camponenes for rol! and pitch aad 
the and Che rate OF chango of the Wenicle's heading are all 
Products Sf fA absgoarte hm The dagtan of atebility tndeor was 
aloo Presented in graphical form in ofder to deresmine ff the 
VWehithe would Pehin under poartive control while manduwiring. 
The moré negative the index, Che fior@ stable tha Wenicle 
metion. In every sec of date the index became more negaiivea 
as rudder aogle wae increased, Therefore. for each ger, the 
atraight liné case wae the beast erable, This samé etfect, 
Ehough not @S pronounced. is presapt for non sero dive plane 
angles. Secommendationa tor fucure research in related areas 
Qneiode tHe evaluation of different contred achemea with 
Thapect to performance, Sone @eamples are the yerification of 
Vehicle peach accuracy, the ability of tha vehielea to fre7ece 
external diatursbancet And remain on courag, 4nd the degree of 
gonfidence in performance apd stability if there are errors 


between the physical and matnemacical modeling of the wehtcte. 
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